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In the present paper, we study finite-temperature phase structure of two- component hard-core 
bosons in a cubic optical lattice. The system that we study in the present paper is an effective model 
for the Bose-Hubbard model with strong on-site repulsions and is called bosonic t-J model. This 
model is a bosonic counterpart of the t-J model for the strongly-correlated electron systems like 
the high-temperature superconducting materials. We study the model by means of path-integral 
methods and Monte- Carlo simulations. We found that this system has a very rich phase structure 
including checkerboard-type "insulating" state, superfluid, phase-separated state, inhomogeneous 
cloudlet state, etc. We are also interested in the possible supersolid phase with both the checker- 
board order and superfluidity and found that additional nearest-neighbor inter-species attractive 
force induces the supersolid state. In the supersolid state, paired superfluid appears in addition 
to the superfluid of single atom. This result gives important insight into mechanism of the high- 
^ ' temperature superconductivity of the cuprate. 

■ 

PACS numbers: 67.85.Hj, 75.10.-b, 03.75.Nt 

I. INTRODUCTION 

In the last decade, system of ultra-cold atoms has been one of the most actively studied area in 
atomic and condensed matter physics. Its various physical properties have been investigated quite 
intensively by both experimental and theoretical methods 111]. As dimensionality and interactions 
between particles are highly controllable and also there are no effects of impurities and defects, 
cold atomic system in an optical lattice is regarded as a "final simulator" for quantum many-body 
systems and it sometimes elevates a purely academic theoretical model to a realistic one. 

It is also expected that obtained knowledge of cold atomic systems gives useful insights into physical 
properties of strongly-correlated many-body systems. In the present paper, we shall investigate a 
system of two-component bosons with strong repulsions. We shall first introduce a bosonic t-J model, 
which is an effective model of the Bose-Hubbard model j2| with strong on-site repulsive interactions 
between atomsfsl-IHI. The original t-J model of fermions is a canonical model for the high-temperature 
(T) superconducting (SC) materials[6]. The fermionic t-J model has been studied quite intensively 
since the discovery of the high- Tc SC, but precise knowledge of its phase structure etc is still lacking 
partially because of difficulty of numerical study of the fermionic system. We expect that study on 
more tractable bosonic counterpart of the t-J model is useful to understand the strongly-correlated 
quantum systems including high- Tc materials. This is one of our motivation to study the bosonic 
t-J model. 

Another motivation of the present study is to see if some exotic state like the supersolid (SS), 
which was theoretically predicted in single-component boson systems, exists in the mult i- component 
system. In 2004, it was reported that the SS state of ^He, i.e., solid with superfiuidity (SF), was 
observed by experiment Soon after this report, possibility of SS in cold atom systems was 
theoretically investigated. Study of the hard-core bosons on a square lattice concluded that the 
realization of the SS was difficult in that system unless a long-range interaction exists [8]. On the 
other hand for the hard-core bosons on a triangular lattice, parameter region for possible SS state has 
been clarified|9|. In the SS of single-component boson, a density wave and SF coexist. In the present 
paper, we shall study the possibility of SS in the two-component hard-core boson systems. It should 
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be remarked that the SS state with paired superfluid and checkerboard order in the bosonic system 
corresponds to the coexisting phase of the antiferromagnetism (AF) and SC in the electron system. 
Recently this AF and SC coexisting phase was really observed for highly homogeneous samples of 
high- Tc materials [lo|. Then it is quite interesting to see how the SS emerges in the present system. 

This paper is organized as follows. In Sec. II, we shall briefly review the derivation of the t-J 
model from the Bose-Hubbard model. Slave-particle operators are introduced to express faithfully 
the hard-core nature of two-component bosons. Then the model is studied by the mean-fleld theory. 
Obtained phase structure at the total fllling factor less than unity is shown. In particular, the phase 
diagram contains the parameter region corresponding to the SS. In Sec. Ill, we study the system by 
means of Monte-Carlo (MC) simulations. The system in grand-canonical ensemble is studied rather 
in detail. At vanishing hopping of bosons, the system reduces to a spin model. As it was already 
shown in the previous paper [5], this system has three phases at flnite T, i.e., antiferromagnetic (AF), 
ferromagnetic (FM) and paramagnetic (PM) phases in terminology of spin model. Then we turn on 
the hopping amplitude to investigate how the SF state appears. All the above three phases result 
in the SF phase as the hopping amplitude is increased. We are particularly interested in phase 
transition from the checkerboard insulating phase (AF phase) to the SF. It is verifled that the phase 
transition is of flrst order with a large hysteresis loop. This result indicates that a coexisting phase 
of the AF solid and the SF may appear. We show that a phase separated (PS) state of the above 
two phases really exists at certain parameters. We investigate whether the Josephson-like tunneling 
of SF through AF solid takes place. 

In Sec. IV, we study the system in canonical ensemble with flxed total particle number. Starting 
from the AF solid, we increase the hopping amplitude and study how phase of the system evolves. 
We found that the system behaves rather diflFerently depending on the density of atoms. At low- 
density region, there exists a single phase transition to the SF. On the other hand for high-density 
region, we found another phase besides those observed in grand-canonical ensemble. This new 
phase is composed of SF cloudlets (flnite magnitude of droplets) in the AF solid background. The 
system evolves as AF crystal cloudlet state PS state SF. Physical properties of each phase 
are investigated by calculating correlation functions and snapshots. In Sec.V, we study eflFects of 
nearest-neighbor (NN) interaction between atoms. In particular we are interested in relation between 
a paired SF (PSF) and the AF order. In the absence of the NN interspecies attractive force, 
the PSF does not appear and the SF of one-body atom has no AF long-range order. However we 
found that by adding the NN interactions between atoms, the state with both the PSF and AF 
order emerges. This result gives an important implication to mechanism of the high- Tc SC of the 
cuprates. Section VI is devoted for conclusion. 

II. MODEL HAMILTONIAN AND SLAVE-PARTICLE REPRESENTATION 

As explained in the introduction, we shall study the t-J model of hard-core bosons in the cubic 
lattice at flnite temperature (r)(5', ^1^. Hamiltonian of the t-J model is derived from the Bose- 
Hubbard model whose Hamiltonian is given as, 

3 3 ^ ^ 

^Hub = - ^ ta{al^-ai + h.c.) - ^ h{bl^^br + h.c.) + U J]](n«r - ^){nbr - ^) 

r,i=l r,i=l r 
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where r denotes site of the cubic lattice, i{= 1,2,3) is the unit vector in the i-th direction (it 
also sometimes denotes the direction index), and and br are boson annihilation operators, ria 
is the number operator of the boson a, and therefore U and Va are inter-species and intra-species 
interactions, respectively. In the present paper, we shall consider the case ta^h <^ 14, V5, and the 
total filling factor of bosons at each site less than unity[13]. Furthermore we set ta = t^ = t and 
Va = V5 and also the chemical potential fica = l^ch- Recently studied ^^Rb-^''Rb atomic system is 
a typical example corresponding to the present model[J^]. More general cases like ta 7^ ^5 will be 
studied in a future publication. 

Effective Hamiltonian in the large on-site repulsion limit, which is called bosonic t-J model, can 
be obtained by the standard methods of expansion in powers oit/U and t/V ^ 



(2.2) 



where pseudo-spin operator 
to the second order of the expansion 



r with Bj. 



ir,brY, a is the Pauh spin matrices, and up 
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(2.3) 



In Eg. ( 12.21) jicr is the chemical potential of hole at site r. In the following discussion, we shall treat 
t, Jz and as free parameters. After obtaining the critical couplings etc, we shall return to the 
expression ( 12.31) and discuss the relation to the Hubbard model. In the system i/tj in Eq.( l2.2l) . the 
physical state at each site r is expanded by three orthogonal basis state vectors {|0), \a) = aj|0), \b) = 
6j|0)}, where |0) is the empty state of the bosons. 
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FIG. 1: Phase structure of model in a— hole density plane obtained by MFT, where a ■ 
flow and SS for supersoild. 



-J±/Jz- S CP stands for supercounter 



Before going into detailed numerical study on the system, it is useful to investigate the groundstate 
properties of the model by the mean-field theory (MFT). We use a variational wave function of bosons 
that has a site-factorized form, 



i*)=n 
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|0). 



(2.4) 



It is straightforward to calculate the expectation value of the Hamiltonian H^j (12.21) in the state |^). 
From Etj = (i^tj), we can obtained a phase diagram of the bosonic t-J model at vanishing T. Details 
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have been already published in the previous paper [5], and therefore here we only show the phase 
diagram in FiglH At the total filling factor unity, the system Hu reduces to the spin model, the 
anisotropic Heisenberg model. We introduce parameter a = —J±/Jz for the anisotropy. At a = — 1, 
phase transition takes place, and for a > —1 antiferromagnetic (AF) phase with Xr=even 7^ Xr=odd 
appears whereas for a < —1 the pure ferromagnetic (FM) phase with Xr = f is realized. For the 
original atomic model, the AF phase corresponds to the state with the checkerboard symmetry of a 
and b atoms. On the other hand, the FM long-range order (LRO) means nonvanishing condensation 
{albr) 0, and sometimes it is called supercounter fiow (SCF) in the original atomic system. As 
holes are doped, SF state with single-atom condensation (a^) = (6^) 7^ takes place. In particular, 
the state with AF+SF is called supersolid (SS) as it has not only the checkerboard symmetry but 
also off-diagonal LRO of the SF. 

For numerical calculation, it is necessary to express the physical state |phys) in the constrained 
Hilbert space faithfully. To this end, we introduce the following slave-particle represent ation[5|, 

ar = (j)lCrl, br = 0jc^2, (2.5) 
(j)l(t)r + cl^Crl + cl2Cr2 " l) |phys) = 0, (2.6) 

where 0^ is a hard-core boson and Cri=i^2 is an ordinary boson. Equation (12.61) is the local constraint 
on the physical state. From Eg. (12.51) . the hard-core boson (j)^ creates the empty state of the original 
bosons ar and 6^, whereas Cri (0^2) creates the state of aj|0) (6j|0)), i.e., 

\0)^4\n), al\0)^cU\Q). bl\0)^cU\n), (2-7) 

where \Q) is the empty state of the slave particles. 

In order to express the local constraint (12.61) in more convenient way, we introduce a CP^ boson 
(the Schwinger boson) z^^, 

Cri = (1 - 0j0r)^r^, (^ = 1, 2) 

(5^i^.^-l)|phys), = 0. (2.8) 

It is easily verified that Eg. (12.61) is satisfied by Eg. (12.81) . The hard-core boson 0^ itself can be 
expressed in terms of another CP-^ boson Wri as follows, 

(l)r = wl^Wri, ^ wl^Wri - l) |phys)^ = 0. (2.9) 

It is easily verified [0^, = 1, and operators (pp^r etc commute with each other. From Eg. (12.91) . 
it is obvious that |0)(^ = wl2\0)w and (/>J|0)(^ = wli\0)w The eflFective Hamiltonian Htj can be easily 
expressed in terms of the CP-^ bosons Zri and Wri. 

III. SYSTEM IN GRAND-CANONICAL ENSEMBLE REVISITED: MONTE-CARLO SIMULATIONS 

In this section, we study the system in the grand-canonical ensemble(GCE). Some results were 
already reported in the previous papersfHt]. In the present section, we shall study the system in an 
inhomogeneous chemical potential jicr and physical properties of a hetero-structure of the AF state 
and superfiuid. In the cold atom system in an optical lattice, the system in the GCE corresponds 
to the case of nearly fiat confining potential. See FigO Generally speaking, around a first-order 
phase transition point at which particle density changes drastically, very subtle adjustment of the 
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chemical potential is required to study the critical behavior of the system. Study of the system in the 
canonical ensemble (CE) with fixed particle density is complementary and helpful for investigation 
of first-order phase transition. In the present paper, the study of the system in the CE will be 
reported in the following section. 

FIG. 2: Physical picture of GCE, in which density of atoms is controhed by chemical potential. 



Partition function of the system in the GCE at temperature T, Zqce, is given by the path-integral 
methods as follows [5], 

Zgce = J [DwDwDzDz\e-^^''^''^''~'^'\ (3.1) 

where /3 = l/{k^T) and the Hamiltonian Htj{w^ z) is obtained from Hu in Eg. ( 12.21) by substi- 
tuting the CP^ variables Zr and Wr for the CP^ operators. We have numerically studied the system 
( 13. Ij) and obtained phase diagram[5[] by calculating the internal energy E and the specific heat (7, 
which are defined by 

E = ^{/3H^j{^,r = 0)), C=^{{f5H,j{^,r = 0)-Ef), (3.2) 

where L is the linear system size. We impose the periodic boundary condition for practical calcula- 
tion, and use the following dimensionless parameters for describing phase diagram etc, 

ci = /3Jz, cs = (3t, a = -J±/Jz. (3.3) 



A. Phase diagram 




FIG. 3: Phase structure of model for t = /Xc = in a — ci plane, where a = —J±/Jz and ci = /SJz. There are three 
phases, PM, AF and FM phases. Physical meaning of each phase is explained in the text. Typical value of particle density 
p = (alar) — (blbr) in each phase is also shown. 



We first consider the case of t = 0, i.e., the system without particle hopping. In Fig|3l we show 
the obtained phase diagram in the a — Ci plane for ficr = and Ci > 0. Details of the calculation 
have been already reported in the previous paper [5,]. There are three phases, paramagnetic (PM), 
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z-component-AF and xy-component-FM phases. At high-T, the system is in the PM phase without 
any long-range order, as it is expected. As T is lowered, phase transition to the AF or FM phase 
takes place. The AF state corresponds to the checkerboard state and the FM state to the state 
with the SCF {a\hr) 7^ 0, as we explained in Sec. II. These phase transitions are of second order. 
In the low-T region, there exists a line of the phase boundary between the AF and FM phases at 
a = — 1 as predicted by the MFT. On the line a = — 1, the pseudo-spin symmetry is enhanced to 
SU(2) as the inter-species and intra-species interactions between atoms are the same, otherwise the 
symmetry is U(l)xZ2, i.e., global rotation of [S^^S^) and reflection —S^. The enhancement 

of the symmetry at a = —1 is seen by the redeflnition of operators as — for r G odd site. 

By this redeflnition, (S^^S^) —{S^^S^) for r G odd site, and then J± —J±- For ferroelectric 
materials, similar transition line to the present a = —1 one is called a morphotropic phase boundary 
and plays a very important role[15]. 
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FIG. 4: Internal energy and hole density as a function of C3 = /3t for ci = 3.0 and a 
to SF is of first order. 
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—0.5. Phase transition from AF state 



Let US turn on the hopping amplitude t. We choose typical three points in the a — Ci plane for each 
of the three phases existing at t = 0, and investigated how the phases change as the value of C3 = /3t 
is increased. We have found that all of the three phases make phase transition to the SF phase. In 
FiglU we show E and hole density for the phase transition from the AF phase to SF. These results 
were obtained by the MC simulations of local-update Metropolis algorithm. More elaborate update 
methods slightly decreases the region of the hysteresis loop[5|]. From the above calculations, it is 
obvious that the phase transition is of flrst order. To verify the flnite superfluidity, we calculated 
the boson correlation function {Bq • 5^), and found it has a nonvanishing correlation as r ^ large 
after the phase transition. We also have measured the spin correlation function and found that the 
SF phase has a FM spin correlation. Physical meaning of the appearance of the FM correlation has 
been explained in the previous paper [5] (see also discussion in Sec. IV). Emergence of the flrst-order 
phase transition stems from the fact that both of the phases separated by the phase transition have 
the LRO's, i.e., the AF and SF+FM orders, respectively. 

Existence of the flrst-order phase transition suggests an interesting possibility of various phases 
in the parameter region of hysteresis loop, e.g.., (i)phase separation (PS) of AF state and SF, (ii) 
supersolid with coexistence of the global AF order and SF, (iii) quantum superposed state of the AF 
and SF cloudlets. We have studied the above problem in both the GCE and the canonical ensembles. 
In the GCE of a constant chemical potential, in which the average number density is not conserved 
as the parameters in the Hamiltonian vary, we have found that in most of the parameter region only 
either of the AF solid or SF is stable in the flat potential with /Hcr = constant, and the values of 
parameters t and /Hcr control which phase appears. In fact, we choose an inhomogeneous state as an 
initial state and update the conflguration by means of the standard Metropolis algorithm. After 10^ 
updates, the whole system tends to become either the AF state with checkerboard order or the SF 
with intermediate homogeneous hole density. Only at very speciflc value of Csc = 9.362, the stable 
PS state of the AF and SF phases appears. This means that the genuine flrst-order phase transition 
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takes place at c^c = 9.362, and the system is composed of immiscible AF and SF domains at the 
phase transition point. 

Here one should remark that the phase transition takes place at ^ 0(1) and /3t ^ 0(10). 

Then ^, y ^ 0{j^). This means that the phase transition observed above in the t-J model also 
exists in the original Hubbard model. 

B. Heterojunction of superfluid and AF crystal 




FIG. 5: Heterojunction of SF and AF state. (Left) Initial state before update. (Right) After 8x10^ sweeps. The central bright 
region is the AF solid state with high particle density, whereas the dark regions represent state of high hole density. Confining 
potential is given by the chemical potential with A = 40 and Ni = 2. 




FIG. 6: Snapshots in a horizontal plane of heterojunction of SF and AF state in Fig|5] (Left) Density of a-atom (Center) Density 
of 6-atom (Right) Density of hole. It is obvious that a atom and b atom reside on each sublattice in the AF region. 

In the previous subsection, we showed that both the AF and SF states can coexist at the first-order 
phase transition region. In the GCE with constant chemical potential, either the AF or SF state is 
realized as a "homogeneous" state except at the genuine critical coupling Csc = 9.362. However in the 
confinement potential of real experiments, the density of the atoms is inhomogeneous. In the central 
region in which the atomic density is higher, the AF state is expected to appear, whereas in the outer 
region the SF state with a lower atomic density is expected to appear. If the inhomogeneous PS 
state is realized around the first-order phase transition point, it resembles to the classical solid-fiuid 
PS state like coexisting state of ice and water. However as the present system is a quantum system, 
one may expect that some interesting and characteristic phenomena of macroscopic quantum system 
occur in such a coexisting state. One example is a kind of the Josephson tunneling in a heterojunction 
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of the AF "crystal" state and SF. In the present case, the both AF state and SF are made of a and 
b atoms. Therefore, to study physical properties of hetero junction of these states is very interesting. 

In this subsection, we study the above problem by means of MC simulations. To specify spatial 
direction of the AF layer, we introduce a "confining potential" in terms of the chemical potential /Licr- 
In the practical calculation, we make an AF layer parallel to the y — z plane at x ^ by adjusting 
chemical potential as ficr = A(> 0) for x = 0,±1,--- ,±A^/, and otherwise zero. As expected, a 
stable AF layer is generated after the MC local update of the system. We show some example in 
Figsl5]and[6l Please notice that we are employing the periodic BC. 




10 20 30 40 50 60 200000 400000 600000 800000 

r i^xkte time 



FIG. 7: Physical properties of heteroj unction of SF and AF state. (Left) Spin correlation (^)- There exists no spin 
correlations in SF region. (Right) Phase difference $12 of Bose condensates in two SF regions as a function of update. It 
fluctuates strongly. 



We first show the spin correlation function for the above heterojunction, 

Gi{r) = {Sr,-Sr,+r), (3.4) 

where the site Tq is located in one of the AF layers and r is a lattice vector perpendicular to the AF 
and SF layers. It is obvious that there exist no spin correlations in the SF regions, whereas in the 
AF regions the checkerboard configurations of a and b atoms are realized. 

In this PS state, it is interesting to investigate if the tunneling of SF current takes place through 
the AF layer. To study this problem, we measured the correlation function of the bosons through 
the intermediate AF layer. If there exists a finite correlation between SF states separated by the AF 
layer, "Josephson tunneling" occurs as in the SC-insulator heterojunction. In FiglTl we show the 
result of the SF correlation $12 that is defined as 

log ( {Bl^- Br,))] ^ (3.5) 



$12 = Im 

where ri (r2) is located inside of the left (right) SF layer in FigO It is obvious that correlation of 
two SF states fiuctuates randomly in the MC updates. We examined some other cases and obtained 
similar results. This result is unexpected because both the AF state and SF are made of a and b 
atoms and then it is expected that atoms move almost freely through the interface. In order words, 
there exists a sharp boundary between the two states and the "classical state", i.e., the AF solid, 
and the quantum state, i.e., SF, have very different properties. 



IV. PHASE SEPARATION AND INHOMOGENEOUS STATES IN CANONICAL ENSEMBLE: 

MONTE-CARLO SIMULATIONS 

In the previous section, we studied the system in the GCE with fixed chemical potential. In 
particular we are interested in the AF phase and SF that are separated by the first-order phase 
transition. In the GCE, in most of the parameter region it is observed that either the homogeneous 
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FIG. 8: Physical picture of CE, in which the average density of atoms is fixed. 



AF state or the homogeneous SF state is stable. However in the previous studies on some related 
systems like hard-core bosons on a triangular lattice, the PS state is observed in the CE even if it is 
not found in the GEC|9|]. This result stems from the fact that the CE corresponds to the physical 
situation of fixed averaged particle density and therefore atoms are assumed to be confined in a 
certain fiat region with a sharp boundary like a box, see FiglHl Then it is interesting to study the 
present system in the CE and see if a SS state or some related state exists in some parameter region. 
To this end, we employ the MC simulations with the standard Metropolis algorithm. The typical 
sweeps for measurement is (30000 ^ 50000) x (10 samples) and the acceptance ratio is 40% ^ 50%. 
Errors are estimated from 10 samples with the jackknife methods. 



By the practical calculation, we found that physical properties of the system, including phase 
structure itself, changes as the average particle density varies. In this subsection, we study the 
system with average atomic density = J^ri^l^r) = Pb = ^ J^ri^l^r) = 15% and the system 
size L = 24. This case corresponds to the low-density region. For small hopping amplitude t and 
a > —1, it is expected that the system at low T tends to separate into the AF solid and "empty" 
region, as the Boltzmann factor dominates the entropy factor at low T. As the hopping amplitude t 
is increased, a phase transition to the SF is expected to takes place as observed by the simulations 
in the GCE in the previous section. This expectation has been actually confirmed by the MC 
simulations. 



We first consider the case Ci = 5 and a = —0.5, relatively large value of J^. For small hopping 
amplitude t, the system exists in the AF phase. As t is increased, a phase transition takes place at 
C3 6.0. See FiglHl i^i which internal energy and specific heat of each term in Hu in (12.21) are shown. 
It is useful to see snapshot of the system before and after the phase transition. See FigHOl It is 
obvious that holes (i.e., atoms) are localized and make some domain for C3 = 0, whereas particles 
are homogeneously distributed for C3 = 8. 

To verify that the observed phase transition is a transition from the AF solid to the SF, we 



A. Low-density region 



1. Large Jz 
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FIG. 9: Internal energy and specific heat of ci and C3 terms as a function of C3 for system with pa — ph — 15% in CE. Results 
show existence of first-order phase transition at C3 ^ 6.0. 
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FIG. 10: Large Jz case. Snapshot of density of holes in the system with pa = pb = 15% before (c3 = 0) and after (c3 = 8) 
the phase transition at cs ^ 6.0. ci = 5.0. Bright (dark) region represents region of high (low) atomic density. In the low- 
hopping region (c3 = 0), system is divided into AF and empty regions. On the other hand in the high-hopping region (c3 = 8), 
homogeneous SF state is realized. The phase transition is of first order. 



calculated the correlation function of the boson operator = (a^, 6r)S 

G'sF(r) = -^5]«-i?.o+.)- (4.1) 
We also calculated the spin correlation function, 

Gsir) = j^J2(Sr,-Sr,+r). (4.2) 

ro 

Results are shown in FigHH It is obvious that the SF phase transition at C3 6.0 accompanies a 
transition from the AF to FM in the spin degrees of freedom. As explained in the previous paper [5|, 
it should be remarked that the FM order in the SF state is a natural outcome of the two-component 
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SF and it should not be regarded as the genuine FM order in the magnetism. In fact, the wave 
function of the a and 6-atomic SF is weU-described by the foUowing MFT-type wave function, 



isF) = n 



1 



COS ■ 



|0), 



(4.3) 



where we have assumed the translational symmetry and recovered the relative phase between a and 
h SF condensates, e"^, which was ignored in Eg. ( 12.41) . By using Eg. ( 14.31) . it is straightforward to 
calculate expectation value of the pseudo-spin operator SV to find that the state |SF) exhibits the 
FM order. 





ci = 5, C3 = ci — 5, C3 

FIG. 11: Correlation functions G^y{t) and Gs{t). 



2. Small Jz 



We also studied the low-density case with Ci = 3 and a = —0.5, relatively small value of J^. 
Calculation of the specific heat C is shown in FigHSl C has a system-size dependence, which 
indicates that the phase transition at C3 c:::^ 5.5 is of second order. We also show snapshots of hole 
density for C3 = 0. Compared with the case of Ci = 5, atoms are distributed rather homogeneously. 
We verified that spin correlation function Gs{r) in that region has no LRO, and therefore the phase 
transition is from the PM to FM+SF. This is the reason why the phase transition turns to second 
order. The result suggests existence of the tricritical point at which the phase transition changes 
from second to first one as Ci is increased. For large J^, there exists a strong attractive force between 
a and b atoms sitting nearest-neighbor sites, and therefore the AF solid tends to be high density. 
On the other hand for small J^, the attractive force is not strong enough to generate the AF solid. 



B. High-density region 

In this subsection, we study the system with average atomic density Pa = Pb = 38% with Ci = 3 
and the system size L = 24. This case is close to that observed by the GCE in the previous section 
and corresponds to the high-density region (see FigEj). Numerical study of the internal energy E 
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FIG. 12: Small Jz case. Specific heat C for pa — ph — 15% with ci = 3 as a function of cs. There exists a second-order phase 
transition at cs 5.5. We also show snapshots for cs = 0. Compared with the case of ci =5, system in the low- hopping region 
is homogeneous. 




FIG. 13: Internal energy and specific heat of each term in i^tj measured in CE as a function of C3. pa — Ph — 38% with ci = 3. 
There are three phase transition s at C3 = 7.6, 8.4, and 14.5. L = 24. 



and the specific heat C indicates the existence of three phase transitions as shown in FigJTSl AU of 
them are of first order as the density of state N{E) has a double-peak shape as shown in Fig JMl 
where N[E) is defined as 

= J dEe-^^N{E). (4.4) 

In Eg. ( 14.41) . J^^ denotes the integration over configurations with a fixed number of atoms (or fixed 
average density of atoms), and N{E) is a function of /3 and the couphng constants. Furthermore 
from 13 to 15^ there is a large hysteresis loop in E indicating existence of a first-order 

phase transition. It is obvious that the system with high particle density has a rather complicated 
phase structure compared with the low-density case. 

Snapshots of the system at various C3 shown in Figll5]are quite useful to understand physical mean- 
ing of the above three phases. At C3 is increased from zero, a transition from the AF checkerboard 
state to a new phase takes place. This new phase is inhomogeneous and composed of immiscible 
AF solid region and SF droplet of mesoscopic magnitude, i.e. SF cloudlet. The appearance of this 
coexisting phase comes from the fact that the AF and SF states tend to have a sharp interface as 
we saw in the previous section. At the second phase transition, a complete PS takes place and the 
system is divided into the AF and SF phases. Width of each phase is determined by the atomic 
density and the hopping amplitude. Finally for large C3, the phase transition into the homogeneous 
SF takes place. 
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FIG. 14: Histograms of N(E). They have double-peak shape, which indicate the existence of the first-order phase transition 
at C3 = 7.6 and 8.4. 




FIG. 15: Snapshot of system at various values of cs with ci = 3. (From left to right)c3 = 7, 8, 10.6 and 16. Bright (dark) region 
corresponds to higher (lower) atomic density. From left to right, AF, cloudlet state, PS and SF. 



It is important to see how the spin correlation function behaves in the above phases. We show 
the results in FigJTGl For C3 = 6, the system exists in the AF phase, and for C3 = 16 the system 
is in the SF phase and the spin correlation has a FM order. In the intermediate regions, the spin 
correlation is a superposition of the AF and FM LRO's. Similar behavior has been observed in other 
systems by the previous study fl6*]. It should be remarked that even in the SP state, the correlation 
function averaged over positions and directions does not exhibit any singular behavior that indicates 
the existence of sharp phase boundary. 

It is interesting to see how SF is realized in each phase. In order to study it, we calculated the 
correlation function of atoms Gsf(^)- The results in Fig JTTI show that in the AF phase the correlation 
tends to vanish very rapidly (no correlation even in the NN pair) , and in the SF phase it has a genuine 
LRO. In the PS state, the correlator seems to have a small but finite LRO, but its finiteness comes 
from the correlation in the SF layers. Then the above results are consistent with the observation by 
the GCE in the previous section, which indicates that the SF correlation terminates at AF layers. 

The second phase existing between C3 ^ 7.6 and 8.4 is a new phase. Counterparts of the other 
three were observed in the previous study of the 2D square lattice at T = o(3]. From the observation 
that the genuine first-order phase transitions separate the second phase from the AF solid state 
and the PS state, it is not a simple mixture of the AF and SF. The most interesting possibility is 
that the phase is a quantum superposition of mesoscopic SF cloudlets in the AF background[17]. We 
verified that the cloudlets change their shape, location and magnitude via MC updates. This fact 
supports the above expectation. Anyway further study is required to clarify physical meaning of 
this interesting phase. 
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FIG. 16: Pseudo-spin correlation for cs = 6, 8, 10 and 16. In the states of mixture of AF solid and SF for cs = 6, 8, 10, 
pseudo-spin correlation is a superposition of AF and FM, whereas in the SF, LRO of FM is realized. 
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FIG. 17: Correlation functions of atoms GsF{r) for cs — 6, 8, 10 and 16. 



V. EFFECTS OF NEAREST-NEIGHBOR INTERACTIONS AND PAIRED SUPERFLUID 



In this section, we shall study effects of NN attractive force between a and b atoms and investigate 
the possibility of the PSF. This problem is closely related to the superconductivity of the cuprates. 
The high- Tc materials are hole-doped AF magnets, and a pair of holes sitting on NN sites feel 
attractive force due to (short-range) AF background. It is widely believed that this attractive force 
is an origin of the superconductivity. The PSF in the bosonic t- J model corresponds to the SC in the 
fermionic model with the condensation of the Cooper pair, whereas the single-atom SF corresponds 
to the metallic state with movable electrons. Then it is very interesting to study if the PSF is 
realized in the present system besides the single-atom SF and to see how the emergent PSF relates 
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to the AF order. 
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FIG. 18: Correlation function of PSF in the cases g = and g = S for pa — ph — 34% and ci = 3. Results indicate existence of 
small but finite LRO for g — 3, whereas no PSF for g — 0. 

Existence of the PSF is examined by measuring the foUowing correlation function, 

GpsF(r) = lJ2{((^lbU^)M^)) - {alao){bUA))- (5-1) 

i 

We show the calculation of Gpsf(^) in the SF state, i.e., Pa = Pb = 34% and Ci = 3, in FigUHl The 
result obviously indicates that there is no finite density of the PSF in the SF of the system Htj. We 
also verified that there exists no PSF in other cases. 

Then we have studied effects of interaction between a pair of atoms on NN sites by adding the 
following term to the Hamiltonian, 

/3^NN = -9Yl {(4^r){bl^^br+^) + (ftj 6.) (a|;+,a,+,)) . (5.2) 

It is obvious that the attractive force i^NN enhances the AF and the PSF. Experimental realization 
of "long-range" interactions like ( 15.21) in cold atom system was recently discussed [il8|]. We have 
investigated the extended t-J model Hu + i^NN with fixed value of g. First for C3 = and positive 
similar phase diagram to that of C3 = 5^ = in Figj3]is obtained but the parameter region of the 
AF state is enlarged as a result of i^NN- 

We are again interested in how the AF state evolves as the hopping parameter C3 is increased. As 
in the case of g = considered in the previous section, we study the low and high-density cases 
separately. 

A. Low-density region 

In this subsection, we consider the case with Pa = Pb = 15% and g = 3. We show E and C as a 
function of C3 in FigJTOl It is obvious that there exists a first-order phase transition at C3 0^ 8.5. For 
small C3, the system is divided into empty state and AF crystal as in the case g = 0^ and after the 
phase transition the homogeneous SF state appears. See FigEOl We calculated Gsf(^), Gpsf(^) and 
Gs{r) in the SF phase, and show the results in Figl2Tl It is obvious that no PSF is generated and 
FM LRO for spin correlation is formed in the SF. This means that effect of the NN attractive force 
is negligibly small in the low-density region. 
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FIG. 19: System of pa — ph — 15%. Internal energy E and specific heat C as a function of cs for g = 3, ci — 3.0 and a ■ 
In C, "specific heat of each term" is shown. 
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FIG. 20: Snapshots cs = (upper panels) and cs = 8 (lower panels) in low-density region with g = 3. For cs = 0, particles a 
and ^ gather in the central region forming the checkerboard pattern, whereas the atoms are distributed rather homogeneously 
for C3 — 8. 



B. High-density region 

Then we studied the system with atomic density Pa = Pb = 34%, i.e., in high-density region. In 
Figl22l we show the specific heat C as a function of C3 for g = 3. The result shows that there exist 
three phase transitions at C3 10.8, 12.5 and 14.5, respectively. We also show the SF and the PSF 
correlation functions in Figl23] (see also FigJTSl). Contrary to the case g = 0^ all SF correlations have 
finite LRO after the first phase transition at C3 0^ 10.8. Small but finite LRO of the PSF exists in 
all phases except the AF phase. 

One may think that this diflFerence in the SF behavior between in the g = and g = 3 cases, in 
particular (non) existence of the PSF, is related to different behavior of spins in the two cases. Then 
we studied the spin correlation and obtained the results shown in FiglSH where 

G^y(^) = ^3j2J2KSlo+r)- (5-3) 

ro i=x,y 

As expected, in the model g = 3 there exists AF correlation in z-component, besides the FM 
correlation in x?/-component. We also studied the snapshots of each phase, which are quite helpful 



17 




FIG. 21: GsF{r) and G'psf(^) (Left panel). Gs{r) (Right panel). In the SF with low-atomic density, no PSF is generated and 
FM RLO for pseudo-spin correlation is formed. 




10 11 12 13 14 15 



C3 



FIG. 22: Specific heat of each term as a function of cs for g = 3, ci — 3.0 and a — —0.5. There are three phase transitions at 
C3 ^ 10.8, 12.5 and 14.5. 



for understanding properties of phases. From these results, we found that the PS takes place even 
in the AF state for C3 < 10.8. However in the present case, the phase is separated into empty phase 
with very low particle density and the AF phase of checkerboard symmetry. In the intermediate 
two phases, the AF solid and SF coexist but they are immiscible. On the other hand in the SF for 
C3 > 14.5, the global domain structure disappears, and a homogeneous state is realized. It is very 
interesting to see if the SS state is realized due to the NN coupling i/^N in Eg. (15.21) . Some previous 
study on the single-band Hubbard like model with a long-range interaction suggested the existence 
of the SS in the square lattice. The above results of the SFs and the spin correlation functions in 
the present case show that the state for C3 > 14.5 is a SS and it appears as a result of the interaction 
i^NN- Snapshots of the SS state in Figl25lshow that atoms are distributed rather homogeneously in 
the most spatial region but there also exist small domains of checkerboard symmetry. The SS state 
is composed of quantum superposition of these inhomogeneous but structureless configurations. It 
is easily seen that the (groundstate) wave function describing the observed SS in the MF level, ^^ss, 
is given as follows, 

l^ss) = n iMbl+i + B,al + B^bl + C] 

rGeven site 

X Yl [Aalbl^. + Boal + Bibl + C]\0), (5.4) 

rGodd site 

where A^B^^Bq and C are constants, and B^. ^ Bo means the existence of the AF order of the 
z-component besides the x^-FM. The wave function (15.41) eloquently shows that the SS state is the 
coherent superposition of the PSF, SF's and holes. Each snapshot shows one component of the 
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FIG. 23: Correlation functions of atoms GsF{r) and G'psF(r) for C3 = 10, 11.5, 13.5 and 15. 
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FIG. 24: Correlation functions of pseudo-spins in each phase of model g — S. There exists AF correlation in ^-component, and 
also FM correlation in x^z-component. 



superposed state ^^ss- 

We expect that the obtained results in this section are quite instructive for the mechanism of 
the high- Tc SC of the cuprates. It is expected that the single-particle SF in the bosonic t-J model 
corresponds to the coherent hopping of electrons in the fermionic t-J model, whereas the PSF 
corresponds to condensation of spin-singlet electron pair at NN sites. In other words, the single- 
particle SF without the PSF corresponds to the (anomalous) metallic phase of the high- Tc materials, 
whereas the single-particle SF with the PSF to the SC state. Effective attractive force appears for a 
pair of holes sitting on NN sites in the AF background, as the NN hole pair break eleven AF bonds, 
whereas two holes separating more than one lattice spacing break twelve AF bonds fio']. It is expected 
that this attractive force is an origin of the high- Tc SC. However the results of the present study 
show that the PSF does not exist in the simple bosonic t-J model, though the phase transition from 
the insulating AF state to the single-particle SF takes place as the hopping amplitude is increased. 
By adding the attractive force for particles sitting NN sites i^NN (15.21) . the PSF is actually realized. 
Furthermore the resultant PSF accompanies the AF order. For the high-Tc cuprates, coexistence 
of the AF and SC was actually observed in the recent clean materials^. This means that the 
attractive force coming from the magnetic terms in i/tj is not strong enough to generate the PSF in 
the present bosonic t-J model. This conclusion may sound rather negative for the idea of AF force 



19 




FIG. 25: Snapshots of density of a-atom (left), 6-atom (center) and hole (right) for C3 = 15 and g = 3. The SS is realized in 
this parameter region. Atoms are distributed rather homogeneously but there exist small regions of the checkerboard pattern 
of the AF. 



for the high- Tc SC. However, we found that the strong enough NN attractive force generates both 
the AF order and the PSF. This indicates a strong correlation between the above two phenomena. 
Anyway, more precise study on difference between effect of coherent fermion hopping and BEC is 
desired. 



VI. CONCLUSION 



In this paper we have studied the three-dimensional bosonic t-J model at finite T. This model 
corresponds to the large on-site repulsion limit of the two-band Bose-Hubbard model. We are par- 
ticularly interested in the case of the fractional filling of atoms and investigated the phase structure. 
We studied the model both in the GCE and CE and found that the model has a very rich phase 
structure. Besides the AF crystal, SF and SP state of the hetero-structure, there exist the inho- 
mogeneous state with SF cloudlet in the AF solid. We found that hetero-structure of the AF solid 
and SF has a sharp phase boundary. From this observation, we give an intuitive picture of the 
hetero-structure that is expected to be observed by experiment in Fig]26l 
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FIG. 26: Intuitive picture of hetero-structure of AF crystal and SF, which is expected to be observed by experiment. 



We also studied the effect of the NN attractive force between a and b atoms. We found that in 
certain parameter region there appears the SS and also PSF. We think that this result gives an 
important insight into mechanism of high- Tc SC of the cuprates. 

In experiment, density and/or mass of two atoms can be varied. Then study of the cases like 
ta 7^ h and/or Pa 7^ Pb is important. This problem is under study and results will be reported in a 
future publication. 



20 



Acknowledgments 

This work was partially supported by Grant-in- Aid for Scientific Research from Japan Society for 
the Promotion of Science under Grant No. 20540264 and No23540301. 



[1] For review, see, e.g., I. Bloch, J. Dalibard, and W. Zwerger, Rev. Mod. Phys.80, 885 (2008); 

M. Lewenstein, A. Sanpera, V. Ahufinger, B. Damski, A. S. De, and U. Sen, Adv. Phy. 56, 243 (2008); 

V.I. Yukalov, Laser Physics 19, 1 (2009). 
[2] D. Jaksch, C. Bruder, J. I. Cirac, C. W. Gardiner, and P. Zoller, Phys. Rev. Lett. 81, 3108 (1998). 
[3] M. Boninsegni, Phys. Rev. Lett. 87, 087201 (2001); Phys. Rev. B 65, 134403 (2002). 
[4] M. Boninsegni and N. V. Prokof'ev, Phys. Rev. B 77, 092502 (2008). 

[5] Y. Nakano, T. Ishima, N. Kobayashi, K. Sakakibara, L Ichinose, and T. Matsui, Phys. Rev. B 83, 235116 (2011); 

Y. Nakano, T. Ishima, N. Kobayashi, T. Yamamoto, I. Ichinose, and T. Matsui, Phys. Rev. A 85, 023617 (2002). 
[6] J.E. Hirsch, Phys. Rev. Lett. 54, 1317 (1985); 

P. W. Anderson, Science 235, 1196 (1987); 

F. C. Zhang and T. M. Rice, Phys. Rev. B37, 3759 (1988). 
[7] E. Kim and M.H.W. Chan, Science 305, 1941 (2004); Nature 427, 6971 (2004). 

[8] P. Sengupta, L. P. Pryadko, F. Alet, M. Troyer, and G. Schmid, Phys. Rev. Lett. 94, 207202 (2005); 

B. Capogrosso, and C. Trefzger, M. Lewenstein, P. Zoher, and G. Pupillo, Phys. Rev. Lett. 104, 125301 (2010). 
[9] S. Wessel and M. Troyer, Phys. Rev. Lett. 95, 127205 (2005). 

[10] H. Mukuda, M. Abe, Y. Araki, Y. Kitaoka, K. Tokiwa, T. Watanabe, A. lyo, H. Kito, and Y. Tanaka, Phys. Rev. Lett. 
96, 087001 (2006). 

[11] A. Kuklov, N. Prokof'ev, and B. Svistunov, Phys. Rev. Lett. 92, 030403 (2004). 

C. Menotti and S. Stringari, Phys. Rev. A 81, 045604 (2010). 

[12] For the bosonic t-J model on a square lattice at T = 0, see Ref.H|. Some comments on three-dimensional system are also 
given there. 

[13] The Bose-Hubbard model at commensurate fillings was studied by, e.g., E. Altman, W. Hofstetter, E. Demler, and M. D. 
Lukin, New J. Phys. 5, 113 (2003); 

S. G. Soyler, B. Capogrosso-Sansone, N. V. Prokof'ev, and B. V. Svistunov, New J. Phys. 11, 073036 (2009). 
[14] S. B. Papp, J. M. Pino, and C. E. Wieman, Phys. Rev. Lett. 101, 040402 (2008). 
[15] Y. Ishibashi and M. Iwata, Jpn. J. Appl. Phys. 37, L985 (1998); 

H. Fu and R. E. Cohen, Nature 403, 281 (2000). 
[16] See for example, K.Aoki, K. Sakakibara, I. Ichinose, and T. Matsui, Phys. Rev. B 80, 144510 (2009). 

[17] For macroscopic quantum superposition of Bose-Einstein condensates, see for example, J. I. Cirac, M. Lewenstein, K. 

M0lmer, and P. Zoller, Phys. Rev. A 57, 1208 (1998). 
[18] K. Coral, L. Santos, and M. Lewenstein, Phys. Rev. Lett. 88, 170406 (2002); 

H.P. Biicher and G. Blatter, Phys. Rev. Lett. 91, 130404 (2003); 

J. M. Sage, S. Sainis, T. Bergeman, and D. DeMille, Phys. Rev. Lett. 94, 203001 (2005). 
[19] I. Ichinose and T. Matsui, Phy. Rev. B 45, 9976 (1992). 



